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We investigate the efTect of white-noise perturbations on chaotic trajectories in open bilhards. We 
focus on the temporal decay of the survival probability for generic mixed-phase-space billiards. The 
survival probability has a total of five different decay regimes that prevail for difi'erent intermediate 
times. We combine new calculations and recent results on noise perturbed Hamiltonian systems 
to characterize the origin of these regimes, and to compute how the parameters scale with noise 
intensity and billiard openness. Numerical simulations in the annular billiard support and illustrate 
our results. 
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In the last decade, experiments and applications 
strengthened the need of considering openness in 
otherwise closed billiards. The relevant dynam- 
ics in such systems is the transient emptying of 
the billiard, instead of time asymptotic properties 
of closed systems. For instance, instead of the 
positive Lyapunov exponents of typical trajecto- 
ries, in open systems the signature of strong chaos 
is an exponential decay of the survival probabil- 
ity inside the billiard. In this paper we perform 
a further step in this asymptotic-to-transient direc- 
tion: we put forward the idea that pre-asymptotic 
decay regimes of the survival probability are well 
defined and physically relevant. We characterize 
the different regimes that dominate the interme- 
diate time dynamics of noise-perturbed billiards, 
and we relate those regimes to the invariant struc- 
tures of the closed phase space. 



I. 



INTRODUCTION 



The idea of considering a hole or leak in an otherwise 
closed billiard is very appealing. This was the physical 
picture used by Pianigiani and Yorke in their influential 
work on conditionally invariant measures [T]. Nowadays 
it is an essential element in the theoretical description of 
experiments in atomic, acoustic, microwave, and optical 
cavities gHl]. 



Another important ingredient considered in this pa- 
per is the effect of random independent and identically 
distributed perturbations (noise) on trajectories. There 
are many good reasons to consider it: test the structural 
stability of the results, facilitate calculations (e.g., ran- 
dom phase approximations), and mimic the effect of well- 
defined physical processes that are too complicated or 
high-dimensional to be modeled in detail (e.g., molecular 
diffusion) . Not surprisingly, the effect of noise is a clas- 
sical problem in dynamical systems. In low-dimensional 
chaotic Hamiltonian systems more than three decades of 
research considered the effect of noise on diffusion and 
anomalous transport [5H7], on the trapping of trajecto- 
ries [HHH], on scattering [T31[T3], etc. 

In this paper we consider the effect of white noise on 
the dynamics of open billiards. We are interested in the 
generic case of billiards with mixed phase space in which 
regions of regular and chaotic motion coexist. We show 
how the survival probability of trajectories inside the bil- 
liard is modified due to noise, and how the transition 
times and parameters of this modified survival probabil- 
ity scale with noise intensity and leak size. The paper has 
a straightforward organization: we start with the closed 
billiard (Sec. [Til), which we leak (Sec. Ill), and perturb 
by noise (Sec. IV). 



II. 



CLOSED BILLIARD 



A. Definition of the dynamics 



* Author to whom correspondence should be sent. E-mail address: 
edugalt@pks.mpg.de 



The annular billiard is defined by two eccentric circles 
with ratio < r < 1 between the two radii and dis- 
tance < (5 < 1 — r between the two centers. Since 
its introduction by Saito et al. three decades ago |15) . 
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FIG. 1. (Color online) The annular billiard with r — 0.3 and S — 0.65. (a) Configuration space showing selected periodic orbits, 
(b) Phase space obtained as a Poincare surface of section at the outer circle. Shown MUPOs: (2, 1) (green triangles), (7, 1) (red 
circles) and KAM island around the period 6 orbit (blue squares). Black dots in (b) correspond to a single chaotic trajectory. 
Three positions of leaks are indicated: loc.i = 0.161, a;c,2 ~ 0.5, and ijJc,3 ~ 0.55 (the shaded region shows the leak I^i^ 3 used 
in Figs. [2] and |3|. 



the annular billiard has been used to investigate differ- 
ent physical phenomena [l6l - fT8] . Here we consider the 
case r = 0.3 and S = 0.65 shown in Fig. llTa). These 
parameters were chosen because numerical simulations 
shown in Fig. [TJb) strongly indicate a mixed phase space 
with the coexistence of regions of regular and chaotic mo- 
tion. Below we discuss the main properties of this partic- 
ular billiard but we emphasize the generality of our main 
results to the generic case of mixed-phase-space billiards. 

The dynamics in the annular billiard can be easily con- 
structed by noting that after colliding with the outer cir- 
cle, trajectories can either collide directly with the outer 
circle or instead first collide once with the inner circle. 
Trajectories that collide with the inner circle fulfill the 
following condition [TS] 



I sin a + (5sin(a — 27rix))| < r. 



(1) 



where a is the angle between the velocity and the nor- 
mal vector at the collision point and w is the normalized 
position of collision at the outer circle (see Fig.flk). The 
significance of the hitting condition in Eq. (IT]) is that it 
can be used to obtain a mapping between two successive 
collisions at the outer boundary [IS] . We investigate this 
billiard map, which corresponds to a Poincare surface of 
section at the outer boundary. For convenience, time 
is counted discretely between successive collisions of the 
map. The collision time (tcoi, assuming constant veloc- 
ity u = 1) in the chaotic component of the annular billiard 
is bounded by min{l — r — S, 2-^1 — {r + Sy} < icoii < 2, 
so that the scalings and the existence of regimes of decay 



remain unaffected by this simplificationM 



B. Phase space components 

The phase space in Birkhoff coordinates (cj, sin a) of 
the closed annular billiard is shown in Fig. fUb) . It can 
be divided in four invariant components (i)-(iv), which 
are built by the trajectories that: 

(i) do not cross the inner circle of radius r + S (sin a > 
r + S) and therefore never satisfy Eq. ([T]). Graph- 
ically they correspond to orbits close to the outer 
boundary of the circle and are called vifhispering 
gallery. For the billiard used in this paper, the 
whispering gallery exists for | sinaj < 0.95, beyond 
the dotted line in Fig. [ij 



(ii) 



cross the circle of radius r + S (sin a < r + 6) but 
never satisfy Eq. (fTl) . These conditions are satisfied 
by periodic orbits that build one-parameter families 
of marginally unstable orbits (MUPOs) [51]. One 
trivial example is the diameter (period p — 2 and 
winding number q — 1 orbit) highlighted in Fig. nl 
In opposite to the whispering gallery, these orbits 
are usually embedded in a chaotic component (see 



^ In ergodic closed systems the correspondence between discrete 
and continuous time is given by the mean collision time (tcoll) = 
^— , where A and d are the area and perimeter of the billiard. 
In open systems this does not hold [191120J. The specific quanti- 
tative results change differently in different parts of the survival 
probability, however, the scalings remain unchanged. 



below) and affect the dynamics of chaotic trajecto- 
ries despite having zero measure, see Ref. [22] for 
a detailed investigation in the annular billiard and 
Refs. [m [^5^ES] in other systems. For the billiard 
used in this paper the following MUPOs {p, q) exist: 
(2, 1), (6, 1), (7, 1), (8, 1), (9, 1), (19, 2), (29, 3), (39, 4), 
(49, 5), (59, 6), (69, 7), (79, 8), (89, 9),.... 

(iii) cross the circle of radius r + S (sin a < r + S), sat- 
isfy Eq. (fTl), but remain close to stable a periodic or- 
bit. These orbits can be periodic, quasi-periodic, or 
even chaotic (confined inside the last quasi-periodic 
circle) and build the so-called Kolmogorov-Arnold- 
Moser (KAM) islands. For the billiard used in this 
paper the most prominent examples are the tra- 
jectories around the period 6 orbit shown as I in 
Fig.g 



(iv) 



cross the circle of radius r + S (sin a < r + S), sat- 
isfy Eq. (fTl), are chaotic, and fill a large component 
of the phase space. Despite the mathematical diffi- 
culties to provide rigorous proofs (see Refs. [27^ |28j 
for rigorous results in particular cases) , it is largely 
believed that the annular billiard with S ^ con- 
tains a large chaotic component, in which a single 
trajectory visits a positive area of the phase space 
[l5Vll8j . This large component is called chaotic 
sea. In Fig. Il|b) this region corresponds to the 
large dotted component. 



III. 



OPEN BILLIARD 



A. Definition of tiie dynamics 

We open the annular billiard by considering a region / 
at the border of the billiard through which the trajec- 
tories escape [H [HI [25l [MSI M IM]. Formally, the 
dynamics of the system with leak M is defined based on 
the dynamics of the closed billiard M as: 



M = 



escape 

M 



for X £ I, 
for X ^ I, 



(2) 



where x— (a;,sina). Notice that, by convention, escape 
occurs only one time step after trajectories hit the leak / 
so that M is defined in /. We are interested in finite but 
small / such that a non-trivial dynamics still exists in the 
billiard. 

Here we consider leaks placed inside the large chaotic 
sea, item (iv) of Sec. II B and we are interested on how 



the trajectories escape from this ergodic component of 
the closed billiard's phase space. The invariant set of 
the open system is the so-called chaotic saddle [33], the 
set of points that never leave the system neither in for- 
ward nor in backward times. In fully chaotic systems the 
chaotic saddle is a zero measure fractal set. In the case 
of mixed phase space system discussed here, the chaotic 



saddle relevant to the escaping trajectories also contains 
a similar hyperbolic component [3TJ [36] but additionally 
it includes a non-hyperbolic component composed by the 



borders of the whispering gallery [region (i) of Sec. II B 
the MUPOs [region (u) of Sec. IIIB], and of the KAM 



islands [region (in) of Sec. II B 



In our simulations of the annular billiard we consider 
leaks / = [wc — Aw, tOc + Aw] x [—A sin a, -l-A sin a] cen- 
tered at three different positions Wc (see Fig. [I]), varying 
Aw, and a fixed A sin a = 2/3. Physically, this configu- 
ration corresponds to a dielectric billiard with refraction 
index 77 — (Asina)"^ = 1.5 (glass) with a perfect mirror 
boundary everywhere except inside the leak, where tra- 
jectories escape for collisions below the critical angle Qc 
with sinofc = l/?7. 



B. Survival probability 

We compute the survival probability P{t) inside the 
billiard by starting an ensemble of trajectories dis- 
tributed according to an initial density po{x). In our 
simulations we consider po{x) to be uniform inside the 
leak of the billiard and elsewhere (for the first iteration 
the closed billiard M is used). Physically, these initial 
conditions correspond to throwing trajectories inside the 
billiard through the leak. Another motivation for using 
this particular po{x) comes from the fact that P{t) in this 
case corresponds exactly to the distribution of Poincare 
recurrence times j31j . The main decay regimes of P{t) re- 
main unaffected by this choice of po{x), but values of the 
exponents and transition times may change [see bullet 
items after Eq. (|3|]. 

Figure 2 shows the decay of the survival probability 
for the particular leak shown in Fig. [Tj We can identify 

mm 



three different regimes of decay I5|)| I 
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deterministic 



(t) 



irregular 


for t < i/. 


ae--'' 


for ti <t <tii, 


bt-P 


for t > tjj. 



(3) 



The dynamics of a typical trajectory escaping in each of 
the regimes in Eq. ([s]) can be associated to the phase 
space structures as: 

9 t < tj irregular: trajectories that collide only a few 
times in the chaotic sea of the closed billiard [region 
(iv) of Sec. II B , the spatial density of survival tra- 
jectories has not converged yet. The exact shape 
in this regime is extremely sensitive to po{x). 

• tj < t < til exponential e~^^: trajectories explore 
the hyperbolic component of the chaotic saddle but 
escape before coming close to the border of the non- 

. Anal- 
kl, the 



hyperbolic components [(i,ii,iii) of Sec. II B 

ogous to the case of fully chaotic systems 

same exponent 7 is observed for different po (2?) [31j . 

• t > til power-law t~^: trajectories get stuck 
close to the non-hyperbolic components of the sad- 
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FIG. 2. (Color online) The survival probability P{t) for the open annular billiard, (a) Logarithmic scale in the y axis; Inset: 
X, y axes in linear scale, magnification for short times, (b) Logarithmic scale in both x, y axes. All regimes and transition times 
in Eq. (|3| are depicted: the dashed line corresponds to the fitting of an exponential with a = 0.72 and 7 — 0.016; and the 
dotted line corresponds to the fitting of a power law with b — 2.492 and /? — 1.65. The transition times are ti = 14 (visual 
inspection) and tu — 1,020. The leak is introduced in the position ljc,3 with ii{I) = 0.01, see Fig. IT] 



die (classical references on such stickiness phenom- 
ena are Refs. [37ti39) ). The asymptotic exponent 
changes from /3 to /?' = /3 — 1 for po{x) nonzero 
at the boundary of the non-hyperbolic components 
(e.g., if po{x) is taken according to the Liouville 
measure restricted to the chaotic sea of the closed 
billiard) [351110]. 

The description above applies to typical trajectories es- 
caping in the corresponding regimes. It is instructive to 
think that for intermediate times, tp < t < tjj, both 
exponential and power-law regimes coexist [53J |3I1 [Ml 



an estimation of the exponential decay as 



P(i)=ae-''*-f 6i-'3, 



(4) 



where i^ is the time needed to approach the non- 
hyperbolic component of the saddle [3T] . 



C. Dependence of the parameters of P{t) on the 
leak 

The exponential decay 7 in Eq. ((sl) can be con- 
sidered a signature of the chaoticity of the map. Fol- 
lowing this reasoning, for small leaks we can approxi- 
mate the escape at each time step by the area of the 
leak relative to the area of the chaotic sea /i(/) := 
Area(/)/Area(chaotic sea). For the billiard considered 
here we found numerically that Area(chaotic sea) « 
0.993 X {r + S) = 0.943. Using phase-space areas cor- 
respond to using the Liouville measure d/i = dud sin a of 
the closed system to approximate properties of the open 
system, and can be shown to be valid for almost all leak 
positions in strongly chaotic systems P^. This leads to 



7 



ln(l - /i(/)) w /i(/) for fi{I) -^ 0. (5) 



Violations of this approximation in the fully chaotic case 
have been extensively discussed in the recent years and 
are particularly large for leaks containing low-period pe- 
riodic orbits of the closed system |5^-^5] . Here we extend 
these previous results and verify the effectiveness of the 
approximation in Eq. ([5| for the intermediate-time ex- 
ponential decay in Eq. H . In Fig. p[a) we compare the 
numerically obtained values to the prediction for differ- 
ent leak sizes centered at three different positions: in the 
chaotic region, around an unstable periodic orbit, and 
around a family of MUPOs. In all cases 7 — > 7* is ob- 
served in the limit of small leaks fi{I) -^ 0, in agreement 
with relation ((sl). For large leak sizes Fig. fig. 3(a) shows 
different deviations of this relation, in agreement with 
the results observed for hyperbohc systems [29H32]| . 

The importance of the value of 7 is that the same value 
is obtained for a broad class of smooth initial densities 
Pq. In fully chaotic systems the requirement is that po 
intersects the stable manifold of the chaotic saddle. Anal- 
ogously, the requirement here is that it intersects the sta- 
ble manifold of the hyperbolic component of the chaotic 
saddles. 

The power-law exponent /3 depends on the proper- 
ties (of the boundary) of the non-hyperbolic sets embed- 
ded in the chaotic component of the phase space (com- 
ponents (i), (ii) and (iii) in the list of Sec. In]). If no KAM 
islands are present, an exponent (3 — 2 can be obtained 
for MUPOs [251I2S1 [34]. The question of whether the 
asymptotic regime has a well-defined and universal power 
law in the generic KAM case is still under investigation 
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FIG. 3. (Color online) Scalings of the parameters of P{t) in Eq. pi with leak size fj,{I) for three different leak locations, (a) 
Intermediate times escape rate 7; The dashed line corresponds to 7* in Eq. (|5|. Inset: dependence of the coefficient a with 
/i(/). (b) Transition time tji, obtained using Eq. ([6|; the dashed line corresponds to a scaling l/fj,{I). The leaks have a fixed 
height A sin a = 2/3 and varying width Alj — >■ 0, which leads to pi(/) — >■ 0. The centers of the leaks are: (jJc,i — 0.161 (in 
a MUPO, black circles), luc.2 = 0.5 (in an unstable periodic orbit, red squares), and uics = 0.55 (in the chaotic sea, green 
diamonds) . 



for the case of area- preserving maps (see Refs. |4T1|42] for 
the latest results that indicate /3 « 1.57). For simplicity 
we write the asymptotic decay as t~^ , but it is meant to 
describe the power-law like behavior usually observed in 
mixed-phase-space systems [37H39] . 



D. Dependence of the transition times of P{t) on 
the leak 



the leak size fi{I) as [3T] 



tii-- 

7 



^^{Iy 



(7) 



for which additional logarithmic corrections apply j23) . 
The scaling in Eq. ([7]) has been confirmed in our numer- 
ical simulations for the three different leak positions, see 
Fig-ib). 



Transition time ti indicates the starting of the expo- 
nential decay. It can be interpreted as a convergence 
time which is proportional to 1/|A'|, where A' is the neg- 
ative Lyapunov exponent of the saddle (the time to relax 
to the hyperbolic component of the saddle along its sta- 
ble manifold). Numerical observations usually show an 
abrupt approach, i.e., the exponential decay provides a 
good description of P{t) after a finite (short) time }31j . 

The transition time tu is defined from Eq. (B as the 
time for which the exponential and power-law contribu- 
tions are equal [31] 



-1*" = btjf ^ p{tii) = 2'^ae 



-/tu 



(6) 



The ratio a/b can be interpreted as the proportion be- 
tween the number of trajectories escaping exponentially 
to the number of trajectories escaping algebraically. It 
depends mainly on the measures of the chaotic and regu- 
lar components of the phase space and therefore it should 
not depend strongly on the measure of the leak /i(/). Un- 
der this assumption we can estimate the scaling of i// on 



IV. 



OPEN NOISY BILLIARD 



A. Definition of the dynamics 

Here we consider additive noise perturbations to the 
trajectories. In the simulations of the annular billiard we 
have implemented at each collision a perturbation to the 
angle a as 



S, 



where S is an independent normal distributed random 
variable with zero mean, (S) = 0, and standard deviation 
a = TT^ (noise strength). In order to prevent the par- 
ticle from leaving the billiard through the border (non- 
physical situation), the noise distribution was truncated 
at a = ±7r/2. Notice that perturbations in the a direc- 
tion are perpendicular to the border of the whispering 
gallery component and to the parameterization of the 
billiard boundary, having therefore a strong impact on 



sliding orbits (a = ±7r/2|^ Based on previous observa- 
tions with different setups [9HT21 [T4] , and in the gener- 
ahty of the arguments below, we believe our results are 
valid for additive white noise in general (provided the 
perturbation in a is nonzero). It is an interesting open 
question whether (and which) modifications are needed 
for multiplicative and colored noise (see Ref. [S]). 



B. Survival probability 



In this description we neglect the effect of noise on 7, 
which has been investigated for a fully chaotic system in 
Ref. jl4|. T he sticky region around MUPOs [region (ii) 
of Sec. II B does not contribute to the enhanced trapping 



regime (to < t < tm) because MUPOs build a zero mea- 
sure set. This means that if noise is added to a system 
in which MUPOs are the only source of stickiness (e.g., 
the Stadium |33] or Drive-belt billiards [35]), we predict 
that the exponential decay will start immediately after 
to (i.e., till = to)- 



In Sec. IIIB the trapping of trajectories inside the bil- 
liard was connected to invariant structures of the deter- 
ministic phase space. The longer the escape time of tra- 
jectories, the closer they approach these invariant struc- 
tures. This leads to a connection between temporal scales 
of the survival probability and spatial scales in the phase 
space. Noise perturbations affect phase-space scales com- 
parable to ^. Based on these arguments we expect that 
for small ^ the survival probability P{t) is modified for 
long times only: 



C. Dependence of the parameters of P{t) on the 
leak 

For small noise perturbations the parameter Prw in 
Eq. ^ can be related to the scaling of the recur- 
rence time distribution of a one-dimensional random walk 
(with step size ^ ^) and is therefore given by ^9^ |43 



Prw 



1 



for ^ ^ 0. 



(10) 



Pit) 



-'deterministic 10^ ^ "^ ^D-) 



P.. 



for t > t]j, 



(8) 



where Pdctcrministic is given by Eq. ^ and to is the tran- 
sition time. Following Refs. [21 [14], the noise perturbed 
survival probability is given by 



Pn. 



.it) 






for to <t < till, 
for t > till. 



(9) 



Figure |4] shows, for ^ = 3 x 10^'', the different decay 
regimes and transition times of the survival probability 
given by Eqs. (l8])-(l9). The dynamics of a typical trajec- 
tory escaping in each of these regimes can be associated 
to the phase space structures of the deterministic dynam- 
ics as: 

• t < tu deterministic: trajectories escape before the 
noise perturbation is noticed, P{t) coincides with 
the C = case. 

• to < t < till enhanced trapping: trajectories enter 
the region corresponding to regular motion of the 
deterministic dynamics [components (i) and (iii) of 
Sec. IIIB] and perform a one-dimensional random 



walk inside it [SHHllTij. 

t > till asymptotic exponential: 
plored all available phase-space. 



trajectories ex- 



^ In fact, trajectories have a tendency of being repelled from sliding 
orbits because the noise perturbation in a is nonzero, truncated 
in a = it7r/2, and added on each collision. Therefore, we do not 
expect these orbits to dramatically affect P{t) or the relation 
between billiard maps and flows in our case. 



In the derivation of this result in Ref. [9] (see also [14]) the 
initial conditions pq{x) were chosen outside KAM islands 
of the deterministic closed billiard, in agreement with the 
case treated here. See Refs. [T0HT2] for the case of poix) 
inside KAM islands. References [II1[I1] also showed that 
if ^^ terms are included, the random walk is biased. 
The scaling in Eq. ( 10 1 is interrupted for long times 



because of the limited region available for the random 
walk in the KAM islands and whispering gallery. In the 
random-walk model this corresponds to adding a reflect- 
ing boundary condition [S] ^J . The exponent 7^ of the 
asymptotic decay can be obtained considering P{t) to 
be a continuous and smooth function around t = tm. 
Evaluating 



dlogiPjt)) 
dt 



Sit t ^ till 



(11) 



for both terms in Eq. ^ and equating then leads to 
tjjj — 0.5/7^. Below we show that tm ^ 1/^^ 



see 



Eq. (14) and Refs. [SHU]], and therefore we obtain 



7« 



1 

till 



(^ for C ^ 0. 



(12) 



This scaling is confirmed in Fig. Isja). Interestingly, we 
observe that for larger noise intensities 75 experiences a 
crossover to 7^ « 7, i.e., 7^ is bounded by the escape 
rate observed for short times which was related to the 
hyperbolic component of the chaotic saddle in Sec. |ll] 
This observation indicates that once the trajectories leave 
the region corresponding to regular motion of the deter- 
ministic dynamics, the hyperbolic component of the sad- 
dle controls their escape. For small noise the process of 
leaving the regular components is slower and therefore 
the scaling (12 1 dominates 7^. For larger noise, the de- 
terministic exponential escape is slower than the escape 
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FIG. 4. (Color online) Survival probability P{t) for the open annular billiard perturbed by noise with intensity ^ = 3 x 10^^. 
(a) Logarithmic scale in the y axis; (b) Logarithmic scale in both x and y axis. All regimes and transition times in Eqs. ([8|-(|9| 
are depicted in (b): fitting of the asymptotic exponential with d = 9.34 x 10~^ and 75 = 1.06 x 10~^ (dot-dashed line in brown); 
power-law decay with c = 0.014 and I3rw = 0.5 (dotted Une in green). The dashed black line corresponds to the ^ = case (see 
Fig. [2|. The transition times were estimated as to = 887 and tm = 15, 158. The leak is as in Fig. |2](a;c,3 with ^{I) — 0.01). 
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from the islands and therefore 75 « 7 is observed. In 
Ref. [10^ a similar but different scaling of 7^ on ^ was 
numerically obtained for the case of random maps and 
for initial conditions taken inside the KAM islands. 



D. Dependence of the transition times of Pit) on 
the leak 



A theory for tu can be found in Ref. [7] and predicts 
that 



Id ^ \lt^ for C ^ 0, 



(13) 



with A ^ 1 related to the scahng of Cantori close to the 
KAM islands. This scaling is valid for to ^ i// because 



for to ~ tii the deterministic trapping around the KAM 
islands has a Umited contribution. 

The final cut-ofi^ time i/// can be estimated from basic 
properties of diffusion motion. The expected distance L 
traveled by a random walker with step-size ~ ^ grows 
as L ^ ^y/i. The time tjji corresponds to the expected 
time for the ^-perturbed trajectories to travel the (fixed) 
distance corresponding to the diameter of the (largest) 
KAM island. Therefore, we estimate 



tiii^e, forego. 



(14) 



Our numerical simulations for the scaling of the tran- 
sition times are shown in Fig. [5lb) . We see that for times 
comparable to i// only a weak dependence of tu on ^ is 
observed. For times larger than tu our results indicate an 
increase of this dependency, consistent with relation ([13 1 . 



The scaling of i/// in Eq. ( 14 1 is confirmed over a larger 
interval of ^. 



combined these results into a complete survival probabil- 
ity which is obtained introducing Eqs. (l3| and ^ into 
Eq. (Is]). For small noise intensities, it contains five dif- 
ferent regimes with four transition times (i/, tu, tm, to) 
and four physically relevant parameters (7,/3,/3ijvy,7^). 
We have discussed how these quantities depend on the 
intensity of the perturbation ^ and on the position and 
size /x(/) of the leak. 

Apart from extending previous results to billiard sys- 
tems, our paper contains several new findings. First, 
we have shown that the (intermediate-times) escape rate 
scales linearly with leak size and is extremely sensitive 
to the location of the leak. These results were previ- 
ously known for fully chaotic systems p9H32] . Second, 
we have shown that the asymptotic exponential decay 7^ 
depends on the noise intensity ^ as 7^ ~ ^^, with a tran- 
sition towards 7^ = 7 for large S, (Fig. Isk). Altogether, 
our results further emphasize the significance of the in- 
termediate time decay regimes of the survival probability 
in weakly chaotic noise-perturbed billiards. 



V. SUMMARY OF CONCLUSIONS 

White noise perturbations have surprising effects on 
the chaotic dynamics of mixed-phase-space Hamiltonian 
systems [5l-fT4]. Here we have shown how these results 
impact the dynamics of open billiards, in which case 
the openness can be controlled systematically and the 
phase space usually contains MUPOs [3T1 122]. We have 
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